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G.  Zhang  and  M.  Gunzburger,  Error  analysis  of  stochastic  collocation  method  for  parabolic  partial 
differential  equations  with  random  input  data',  SIAM  Journal  on  Numerical  Analysis,  50,  2012,  1922-1940. 

A  stochastic  collocation  method  for  solving  linear  parabolic  partial  differential  equations  with  random 
coefficients,  forcing  terms,  and  initial  conditions  was  analyzed.  The  input  data  were  assumed  to  depend  on 
a  finite  number  of  random  variables.  Unlike  previous  analyses,  a  wider  range  of  situations  were  consid¬ 
ered,  including  input  data  that  depend  nonlinearly  on  the  random  variables  and  random  variables  that  are 
correlated  or  even  unbounded.  We  provided  a  rigorous  convergence  analysis  and  demonstrated  the  expo¬ 
nential  decay  of  the  interpolation  error  in  the  probability  space  for  both  finite  element  semidiscrete  spatial 
discretizations  and  for  finite  element,  CrankNicolson  fully  discrete  space-time  discretizations.  Ingredients 
in  the  convergence  analysis  include  the  proof  of  the  analyticity,  with  respect  to  the  probabilistic  parameters, 
of  the  semidiscrete  and  fully  discrete  approximate  solutions.  A  numerical  example  is  provided  to  illustrate 
the  analyses. 

We  provided  the  first  rigorous  analysis  of  stochastic  collocation  methods  for  parabolic  PDE  with 
random  inputs. 

J.  Ming  and  M.  Gunzburger,  Efficient  numerical  methods  for  stochastic  partial  differential  equations 
through  transformation  to  equations  driven  by  correlated  noise',  International  Journal  for  Uncertainty  Quan¬ 
tification,  3  2013, 321-339. 

A  procedure  was  provided  for  the  efficient  approximation  of  solutions  of  a  broad  class  of  stochastic  par¬ 
tial  differential  equations  (SPDEs),  that  is,  partial  differential  equations  driven  by  additive  white  noise.  The 
first  step  was  to  transform  the  given  SPDE  into  an  equivalent  SPDE  driven  by  a  correlated  random  process, 
specifically,  fhe  Ornsfein-Uhlenbeck  process.  This  allowed  for  fhe  use  of  fruncafed  Karhunen-Eove  expan¬ 
sions  and  sparse-grid  mefhods  for  fhe  efficienl  and  accurafe  approximafion  of  fhe  inpuf  sfochaslic  process 
in  terms  of  few  random  variables.  Defails  of  fhe  procedure  were  given  and  ifs  efficacy  was  demonsfrafed 
fhrough  compufafional  experimenfs  involving  fhe  sfochaslic  heal  equation  and  fhe  sfochaslic  Navier-Slokes 
equations. 

We  provided  a  very  efficient  means  for  treating  PDEs  driven  by  white  noise  which  typically  require  a 
much  larger  number  of  solutions  of  the  PDE  to  build  accurate  statistical  ensembles. 

G.  Zhang,  M.  Gunzburger,  and  W.  Zhao,  A  sparse-grid  method  for  multi-dimensional  backward 
stochastic  differential  equations'.  Journal  of  Compufafional  Malhemalics,  31,  2013,  221-248. 

A  sparse-grid  melhod  for  solving  mulli-dimensional  backward  sfochaslic  differenlial  equations  (BS- 
DEs)  based  on  a  mulli-slep  time  discrelizalion  scheme  previously  developed  by  some  of  fhe  aulhors  was 
developed.  In  fhe  mulli-dimensional  spalial  domain,  i.e.  fhe  Brownian  space,  fhe  condilional  malhemal- 
ical  expeclalions  derived  from  fhe  original  equalion  were  approximated  using  sparse-grid  Gauss-Hermile 
quadralure  rule  and  (adaptive)  hierarchical  sparse-grid  interpolation.  Error  eslimales  were  proved  for  fhe 
proposed  fully-discrele  scheme  for  mulli-dimensional  BSDEs  wilh  cerlain  lypes  of  simplified  generator 
funclions.  Einally,  several  numerical  examples  were  provided  lo  illuslrale  fhe  accuracy  and  efficiency  of  our 
scheme. 

In  many  situations,  BSDEs  provide  an  efficient  means  for  solving  deterministic  PDEs  by  solving  a 
simpler,  albeit  stochastic,  differential  equation.  In  this  work,  we  provided  a  novel  and  very  efficient  ( relative 
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to  existing  approaches)  for  solving  BSDEs. 

M.  D’Elia  and  M.  Gunzburger,  Coarse- grid  sampling  interpolatory  methods  for  approximating  Gaus¬ 
sian  random  fields',  SIAM/ASA  Journal  on  Uncertainty  Quantification,  1,  2013,  270-296. 

Random  fields  can  be  approximated  using  grid-based  discretizations  of  their  covariance  functions  fol¬ 
lowed  by,  e.g.,  an  eigendecomposition  (i.e.,  a  KarhunenLo‘eve  expansion)  or  a  Cholesky  factorization  of 
the  resulting  covariance  matrix.  We  considered  Gaussian  random  fields  and  analyzed  the  efficiency  gains 
obtained  by  using  low-rank  approximations  based  on  constructing  a  coarse  grid  covariance  matrix,  followed 
by  either  an  eigendecomposition  or  a  Cholesky  factorization  of  that  matrix,  followed  by  interpolation  from 
the  coarse  grid  onto  the  fine  grid.  The  result  is  coarser  sampling  and  smaller  decomposition  or  factoriza¬ 
tion  problems  than  that  for  full-rank  approximations.  Using  one-dimensional  experiments,  we  examined  the 
relative  merits,  with  respect  to  accuracy  achieved  for  the  same  computational  complexity,  of  the  different 
low-rank  approaches.  We  found  that  interpolation  from  the  coarse  grid  combined  with  the  Cholesky  factor¬ 
ization  of  the  coarse  grid  covariance  matrix  yields  the  most  efficient  approach.  Of  course,  these  approaches 
rely  on  knowing  the  correlation  function  of  the  random  field.  When  that  function  is  not  known  (which  is 
often  the  case),  we  explore  the  use  of  only  interpolation  from  a  coarse  grid  onto  a  fine  grid  as  an  alternative 
to  what  is  usually  done,  which  is  assuming  a  covariance  function.  Computational  tests  performed  on  one- 
and  two-dimensional  uniform  and  nonuniform  grids  indicated  that  this  approach,  which  does  not  require  any 
matrix  factorizations  or  decompositions,  resulted  in  correlations  that  can  be  very  similar  to  those  that  result 
from  assuming  well-known,  e.g.,  exponential-type,  correlation  functions. 

We  developed  grid-based  methods  for  dicretizing  correlated  random  fields  which  can  be  implemented 
in  much  the  same  simple  way  as  for  the  ubiquitously  used  grid-based  methods  for  white  noise  fields. 

A.  Labovsky  and  M.  Gunzburger,  An  efficient  and  accurate  method  for  the  identification  of  the  most 
influential  random  parameters  appearing  in  the  input  data  for  PDFs',  SIAM/ASA  Journal  on  Uncertainty 
Quantification,  2,  2014,  82-105. 

Cut-HDMR  expansions  (also  referred  to  as  anchored- ANOVA  expansions)  have  often  been  used  to 
represent  multivariate  functions  in  high  dimensions  because  they  can  be  used  to  identify  unimportant  vari¬ 
ables.  Past  efforts  in  this  direction  have  examined  only  the  separate  influence  of  each  variable.  However, 
we  showed,  using  simple  examples,  that  variables  that  have  small  separate  influences  can  have  large  inter¬ 
actions,  and  thus  those  variables  should  not  be  ignored.  We  developed  a  methodology  for  determining  the 
importance  of  variables  by  examining  both  their  separate  and  pairwise  effects.  As  a  result,  not  only  were 
all  unimportant  terms  omitted  from  the  cut-HDMR  expansion,  but  all  important  terms  were  retained.  This 
is  in  contrast  to  existing  methods  that  can  omit  important  pairwise  interaction  terms.  The  application  and 
effectiveness  of  the  new  methodology  were  demonstrated  for  a  nonlinear  system  of  partial  differential  equa¬ 
tions  having  random  inputs;  specifically,  we  considered  a  magnetohydrodynamics  setting  which  also  serves 
to  illustrate  that  a  realistic  problem  can  indeed  involve  random  input  parameters  that  have  small  individual 
influences  but  large  pairwise  influences.  In  such  settings,  the  new  method  is  not  only  computationally  at¬ 
tractive,  but  it  is  the  only  method  that  can  correctly  identify  all  important  individual  effects  and  pairwise 
interactions.  Also  considered  is  a  possible  further  application  of  the  new  methodology,  namely  reducing  the 
cost  of  sparse-grid  approximations  of  quantities  of  interest  that  depend  on  the  solution  of  a  partial  differential 
equation. 

Reducing  the  number  of  random  parameters  that  one  must  deal  with  in  problems  described  by  PDFs 
with  random  inputs  is  of  paramount  importance  for  reducing  the  very  high  costs  associated  with  uncertainty 
quantification  for  such  cases.  In  this  work,  we  provided  the  first  means  for  correctly  identifying  this  parame¬ 
ters  that  are  important  not  just  in  their  separate  effect,  but  also  in  their  pairwise  interactions.  An  algorithm 
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to  do  this  was  developed  and  tested  and  applied  to  an  MHD  problem. 

M.  Gunzburger,  C.  Webster,  and  G.  Zhang,  Stochastic  finite  element  methods  for  partial  differen¬ 
tial  equations  with  random  input  data',  Acta  Numerica,  23,  2014,  521-650. 

Acta  Numerica  has  for  over  20  years  been  the  leading  outlet  for  publishing  timely  review  articles  on 
all  aspects  of  computational  mathematics  and  numerical  analysis.  It  has  been  especially  used  to  acquaint 
readers  with  emerging  areas  of  research.  I  was  invited  to  prepare  and  article  on  uncertainty  quantification 
which  is  absolutely  an  emerging  areas  in  mathematics  and  applications.  The  result  was  this  130  page 
article  with  provides  a  wealth  of  rigorous  theoretical  and  practical  useful  knowledge  about  uncertainty 
quantification  involving  PDFs. 

C.  Webster,  G.  Zhang,  and  M.  Gunzburger,  An  adaptive  sparse-grid  iterative  ensemble  Kalman 
filter  approach  for  parameter  field  estimation',  International  Journal  of  Computer  Mathematics,  91,  2014 
798-817. 

A  new  sparse-grid  (SG)  iterative  ensemble  Kalman  filter  (lEnKF)  approach  was  developed  for  estimat¬ 
ing  spatially  varying  parameters.  The  adaptive  high-order  hierarchical  sparse-grid  (aHHSG)  method  was 
adopted  to  discretize  the  unknown  parameter  field.  An  lEnKF  was  used  fo  explore  fhe  paramefer  space 
and  esfimafe  fhe  surpluses  of  fhe  aHHSG  inferpolanf  af  each  SG  level.  Moreover,  fhe  esfimafed  aHHSG 
inferpolanf  on  coarser  levels  is  employed  fo  provide  a  good  initial  guess  of  fhe  lEnKF  solver  for  fhe  approx¬ 
imation  on  fhe  finer  levels.  The  mefhod  is  demonsfrafed  in  esfimafing  permeability  field  in  flows  fhrough 
porous  media. 

By  combining  hierarchical  sparse  grids  with  iterative  ensemble  Kalman  filtering,  both  effective  in 
some  respects  in  their  own  right,  we  have  developed  a  very  efficient  method  for  parameter  field  estimation. 

M.  Gunzburger  and  C.  Webster,  Uncertainty  quantification  in  the  partial  differential  equation  set¬ 
ting',  fo  appear  in  Mafhemalical  Infelligencer. 

The  Mathematical  Intelligencer,  published  by  Springer,  is  the  leading  organ,  on  the  worldwide  stage, 
for  providing  information  about  developments  and  trends  in  all  of  mathematics.  1  was  invited  to  prepare  this 
article  on  uncertainty  quantification,  arguably  the  hottest  trend  in  applied  mathematics  and  statistics. 

G.  Zhang,  C.  Webster,  and  M.  Gunzburger,  A  hyper-spherical  adaptive  method  for  discontinuity 
detection',  submitted. 

We  developed  and  analyzed  a  hyper-spherical  adaptive  hierarchical  sparse-grid  mefhod  for  defecfing 
jump  disconfinuifies  of  funcfions  in  high-dimensional  spaces.  The  mefhod  was  mofivafed  by  fhe  fheorefical 
and  compufafional  inefficiencies  of  well-known  adaptive  sparse-grid  mefhods  for  disconfinuify  defection. 
Our  novel  approach  consfrucfed,  using  a  hyper-spherical  fransformafion,  a  funclion  represenfafion  of  fhe 
disconfinuify  hyper-surface  of  an  A^-dimensional  disconfinuous  quantify  of  inferesf.  Then,  a  sparse-grid 
approximation  of  fhe  fransformed  funcfion,  builf  in  fhe  hyper-spherical  coordinafe  system,  whose  value 
af  each  poinf  was  esfimafed  by  solving  a  one-dimensional  disconfinuify  defection  problem.  Due  fo  fhe 
smoofhness  of  fhe  hyper-surface,  fhe  new  fechnique  could  identify  jump  disconfinuifies  wifh  significanfly 
reduced  compufafional  cosf,  compared  fo  exisfing  mefhods.  Moreover,  hierarchical  accelerafion  fechniques 
were  also  incorporafed  fo  furlher  reduce  fhe  overall  complexity.  Rigorous  complexify  analyses  of  fhe  new 
mefhod  were  provided  as  were  several  numerical  examples  fhaf  illusfrales  fhe  effecliveness  of  fhe  approach. 

The  notoriously  difficult  problem  of  discontinuity  detection  with  respect  to  dependences  on  random 
parameters  is  the  crucial  task  one  faces  in  many  risk  assessment,  policy  decision,  and  other  uncertainty 
quantification  settings.  Or  approach,  which  does  not  make  the  futile  attempts  trying  to  identify  discontinu- 
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ities  in  parameter  space,  but  rather  does  so  in  a  transformed  space,  provides  a  novel  and  more  effective 
approach  for  discontinuity  detection. 

A.  Teckentrup,  P.  Jantsch,  C.  Webster,  and  M.  Gunzburger,  A  multilevel  stochastic  collocation 
method  for  partial  differential  equations  with  random  input  data',  submitted. 

Stochastic  collocation  methods  for  approximating  the  solution  of  PDEs  with  random  input  data  (e.g., 
coefficients  and  forcing  terms)  suffer  from  the  curse  of  dimensionality  whereby  increases  in  the  stochastic 
dimension  cause  an  explosion  of  the  computational  effort.  We  developed  and  analyzed  a  multilevel  version 
of  the  stochastic  collocation  method  that,  as  was  the  case  for  multilevel  Monte  Carlo  (MLMC)  methods, 
used  hierarchies  of  spatial  approximations  to  reduce  the  overall  computational  complexity.  In  addition,  our 
proposed  approach  utilized,  for  approximation  in  stochastic  space,  a  sequence  of  multi-dimensional  inter- 
polants  of  increasing  fidelity  which  were  then  used  for  approximating  statistics  of  the  solution  as  well  as 
for  building  high-order  surrogates  featuring  faster  convergence  rates.  A  rigorous  convergence  and  compu¬ 
tational  cost  analysis  of  the  new  multilevel  stochastic  collocation  method  was  provided,  demonstrating  its 
advantages  compared  to  standard  single-level  stochastic  collocation  approximations  as  well  as  MLMC  meth¬ 
ods.  Numerical  results  were  provided  that  illustrate  the  theory  and  the  effectiveness  of  the  new  multilevel 
method. 

Mutilevel  methods,  e.g.,  multi  grid  methods,  have  historically  been  shown  to  provide  avenues  for 
great  improvement  in  the  efficiency  of  numerical  algorithms.  Our  approach  that  combines  hierarchies  of 
spatial  approximation  and  sparse-grid  approximations  in  parameter  space  greatly  improve  the  efficiency  of 
solving  PDEs  with  random  inputs. 

M.  Gunzburger  and  A.  Teckentrup,  Optimal  point  sets  for  total  degree  polynomial  interpolation  in 
moderate  dimensions',  submitted. 

Lagrange  interpolation  by  total  degree  polynomials  in  moderate  dimensions  was  the  focus  of  this 
project.  In  particular,  we  were  interested  in  characterizing  the  optimal  choice  of  points  for  the  interpola¬ 
tion  problem,  where  we  defined  fhe  opfimal  inferpolafion  poinfs  as  fhose  which  minimize  fhe  Lebesgue 
consfanf.  We  developed  a  novel  algorifhm  for  numerically  compufing  fhe  location  of  fhe  optimal  poinfs, 
which  was  independenf  of  fhe  shape  of  fhe  domain  and  does  nof  require  compufafions  wifh  Vandermonde 
mafrices.  We  performed  a  numerical  sfudy  of  fhe  growfh  of  fhe  minimal  Lebesgue  consfanf  wifh  respecf  fo 
fhe  degree  of  fhe  polynomials  and  fhe  dimension,  and  so  produced  fhe  lowesf  values  known  as  yef  of  fhe 
Lebesgue  consfanf  in  fhe  unif  cube  and  fhe  unif  ball  in  up  fo  10  dimensions. 

For  passive  approximation,  i.e.,  barring  the  use  of  adaptive  strategies,  total  degree  polynomial  in¬ 
terpolation  provides  the  most  efficient  means  for  providing  approximations  of  a  desired  precision.  We  have 
overcome  the  difficulties  faced  with  producing  a  good  set  of  interpolation  points  and,  in  fact,  produced  a 
means  for  determining  optimal  interpolations  points. 

H.-W.  VAN  Wyk,  M.  Gunzburger,  J.  Burkardt,  and  M.  Stoyanov,  Power-law  noises  over  general 
spatial  domains  and  on  non-standard  meshes',  submitted. 

Power-law  noises  abound  in  nafure  and  have  been  observed  exfensively  in  bofh  lime  series  and  spalially 
varying  environmenlal  paramelers.  Allhough  recenl  years  have  seen  fhe  exfension  of  Iradilional  slochas- 
lic  parlial  differential  equalions  fo  include  syslems  driven  by  fraclional  Brownian  motion,  spatially  dis- 
Iribufed  scale-invariance  has  received  comparatively  little  attention,  especially  for  paramelers  defined  over 
non-sfandard  spalial  domains.  Our  work  focused  on  fhe  generalizalion  of  power-law  noises  fo  arbilrary 
multidimensional  spalial  domains  by  eslablishing  Iheir  Iheir  Iheorelical  underpinnings  as  well  as  addressing 
Iheir  numerical  simulation.  Three  compulalional  algorilhms  were  developed  for  efficienlly  generaling  Iheir 
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sample  paths,  accompanied  by  numerous  numerical  illustrations. 

White  noise  random  fields  are  in  ubiquitous  use  in  practice  for  modeling  uncertainty  in  complex 
systems,  despite  the  fact  that  the  have  infinite  variance.  On  the  other  hand,  power-law  noises  have  been 
often  observed  but  are  difficult  to  simulate  in  eve  two  and  three  spatial  dimension.  We  have  developed 
efficient  means  for  doing  so. 

G.  Zhang,  W.  Zhao,  C.  Webster,  and  M.  Gunzburger,  Numerical  solution  of  backward  stochastic 
differential  equations  with  jumps  for  a  class  of  nonlocal  diffusion  problems',  submitted. 

We  developed  a  novel  numerical  approach  for  linear  nonlocal  diffusion  equations  with  integrable  ker¬ 
nels,  based  on  the  relationship  between  the  backward  Kolmogorov  equation  and  a  class  of  backward  stochas¬ 
tic  differential  equations  (BSDEs)  driven  by  Levy  processes  with  jumps.  The  nonlocal  diffusion  problem 
under  consideration  was  converted  into  a  BSDE,  for  which  numerical  schemes  were  developed  and  applied 
directly.  The  most  significant  advantage  of  our  approach  is  that  the  BSDE  can  be  solved  independently 
along  each  trajectory  of  the  underlying  Levy  processes,  and  therefore,  completely  avoids  the  difficulties 
associated  with  solving  sequences  of  linear  systems  required  by  traditional  numerical  approximations.  In 
addition,  the  inherent  independence  of  our  procedure  allows  for  embarrassingly  parallel  implementations 
and  also  enables  adaptive  approximation  techniques  to  be  incorporated  in  a  straightforward  fashion.  Rigor¬ 
ous  error  analysis  of  the  new  method  were  developed  as  were  several  numerical  examples  that  illustrated  the 
effectiveness  and  efficiency  of  fhe  proposed  approach. 

Standard  discretization  of  nonlocal  diffusion  problems  result  in  matrices  that  significantly  less  sparse 
than  the  corresponding  discretizations  of  PDFs.  Our  approach  completely  avoids  the  need  to  solve  matrix 
problems,  so  that  it  provides  and  efficient  means  for  solving  nonlocal  diffusion  problems. 

G.  Zhang,  D.  Lu,  M.  Ye,  M.  Gunzburger,  and  C.  Webster,  An  adaptive  sparse-grid  high-order 
stochastic  collocation  method  for  Bayesian  inference  in  groundwater  reactive  transport  modeling',  Wafer 
Resources  Research,  49,  2013,  1-22. 

Bayesian  analysis  has  become  vifal  fo  uncerfainly  quanlificalion  in  groundwater  modeling,  buf  ifs  ap¬ 
plication  has  been  hindered  by  fhe  compufafional  cosf  associafed  wifh  numerous  model  executions  required 
fo  explore  fhe  posterior  probabilify  densify  function  (PPDE)  of  model  parameters.  This  is  particularly  fhe 
case  when  fhe  PPDE  is  esfimafed  using  Markov  Chain  Monfe  Carlo  (MCMC)  sampling.  We  developed  a 
new  approach  fo  improve  fhe  compufafional  efficiency  of  Bayesian  inference  by  consfrucfing  a  surrogafe  of 
fhe  PPDE,  using  an  adaptive  sparse-grid  high-order  stochastic  collocation  (aSG-hSC)  mefhod.  Unlike  pre¬ 
vious  works  using  firsf-order  hierarchical  bases,  our  approach  utilized  a  compacfly  supporfed  higher-order 
hierarchical  basis  fo  consfrucf  fhe  surrogate  sysfem,  resulting  in  a  significanf  reduction  in  fhe  number  of 
required  model  executions.  In  addition,  using  fhe  hierarchical  surplus  as  an  error  indicator  allowed  locally 
adapfive  refinemenf  of  sparse  grids  in  fhe  parameter  space,  which  furlher  improved  compufafional  efficiency. 
To  efficienfly  build  fhe  surrogate  sysfem  for  fhe  PPDE  wifh  mulfiple  significanf  modes,  opfimizafion  fech- 
niques  were  used  fo  identify  fhe  modes,  for  which  high-probabilify  regions  were  defined  and  componenfs 
of  fhe  aSG-hSC  approximation  were  consfrucfed.  Affer  fhe  surrogafe  was  defermined,  fhe  PPDE  can  be 
evaluated  by  sampling  fhe  surrogate  sysfem  direcfly  wifhouf  model  execufion,  resulting  in  improved  ef¬ 
ficiency  of  fhe  surrogafe-based  MCMC  compared  wifh  conventional  MCMC.  The  developed  mefhod  was 
evaluated  using  fwo  synfhefic  groundwafer  reacfive  fransporf  models.  The  firsl  example  involved  coupled 
linear  reacfions  and  demonsfrafed  fhe  accuracy  of  our  high-order  hierarchical  basis  approach  for  approxi¬ 
mating  high-dimensional  posteriori  disfribufion.  The  second  example  was  highly  nonlinear  because  of  fhe 
reacfions  of  uranium  surface  complexafion,  and  demonsfrafed  how  fhe  iferafive  aSG-hSC  mefhod  was  able 
fo  capfure  multimodal  and  non-Gaussian  fealures  of  PPDE  caused  by  model  nonlinearify.  Bofh  experimenfs 
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showed  that  aSG-hSC  is  an  effective  and  efficient  tool  for  Bayesian  inference. 

Although  our  novel  approach  was  developed  in  the  context  of  groundwater  flow  modeling,  the  algo¬ 
rithmic  innovations  apply  to  many  other  Baysean  inference  setting. 

J.  Choi  and  M.  Gunzburger,  Approximation  and  application  of  the  Musiela  stochastic  PDE  in  forward 
rate  models',  International  Journal  of  Computer  Mathematics,  89,  2012,  1269-1280. 

We  considered  the  Musiela  equation  of  forward  rates;  this  is  a  hyperbolic  stochastic  partial  differential 
equation.  A  weak  formulation  of  the  problem  using  the  streamline  upwind  PetrovGalerkin  method  was 
analysed.  Error  analyses  of  the  method  yielded  estimates  for  the  convergence  rates.  Computational  examples 
were  provided  that  illustrated  not  only  the  discretization  methods  used,  but  the  type  of  results  relevant  to 
bond  pricing  that  can  be  obtained  from  the  equation. 

Our  rigorously  justified  methodologies  provide  a  superior  means  for  approximating  forward  rate 

models 
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